Magnetic influence on classical dispersion 

T. N. C. Mendes*, Reinaldo de Melo e Souza^, C. Farina* 
Instituto de Ffsica - UFRJ - CP 68528 
Rio de Janeiro, RJ, Brasil - 21945-970. 

Abstract 

We discuss the Lorentz model for dispersion and absorption of radiation 
in dilute, linear and isotropic materials. Initially, with the purpose of making 
the paper as self-contained as possible, we reproduce the usual calculations con- 
cerning the interaction between the charged material oscillators and the electric 
field of the incident radiation, obtaining the main behavior of the reactive and 
dissipative electromagnetic properties of the materials. Thereafter, we also in- 
clude the magnetic contribution of the Lorentz force to the equation of motion 
of the oscillators up to first order in v/c, which leads to some interesting results, 
like the approximately linear dependence of the refraction index with the radi- 
ation intensity and the appearance of a second region of anomalous dispersion 
around half the natural frequencies of the material. 
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1 Introduction 



As described by A. Pais, [lj, Kirchoff's 1859 pioneering work about the search of the 
spectral function of thermal radiation raised the discussion of the interaction between 
matter and radiation, an issue that, since then, has shown to be of great importance 
in many branches of physics, like quantum field theory, condensed matter physics and 
optical physics, to mention just a few. In this work, we discuss the Lorentz model, one 
of its most prodigious child. Born in the last years of the 19-th century, this model 
treats matter as posessing stable electronic charge configurations, held together by 
harmonic potentials, and that interact with electromagnetic fields of a given radiation. 
This model achieves amazing results describing optical and dissipative properties of 
material media (for instance, their dielectric constants, magnetic permeabilities and 
absorption/dispersion factors), specially those with low refractive index. Usually, 
when discussing classical dispersion theory, only the electric contribution is considered 
since this contribution is c/v times larger than the magnetic one. These issues are 
often covered by undergraduation textbooks in classical electromagnetism [21 [3] , but 
are discussed in more advanced textbooks as well [I]. 

Here, we will show that the inclusion of the usually ignored magnetic term leads 
to quite interesting results as, for instance, the appearence of a second region of 
anomalous dispersion. We also present some numerical analysis in order to see if 
experimental verifications of the new effects discussed in this paper are plausible, 
particulary for electromagnetic radiation of high enough intensity. 

We finish this brief introduction with a comment to motivate the search for other 
possibilities of ocurrence of anomalous dispersion. Besides their intrinsic importance, 
it has been known for a long time that group velocity in regions of anomalous disper- 
sion close to an absorption line may exceed the speed of light in vacuum and, under 
some circunstances, may become infinite or even negative [5j |6|. For a recent and 
pedagogical discussion on negative group velocity see [7] and references therein. 

2 Lorentz's model 

For future convenience, and to make this article as self-contained as possible, we start 
with a brief review of the Lorentz model for the electric permittivity e(u) of linear 
and isotropic media. The interaction between charges and fields in this model is given 
by the Lorentz force. For the j-th charged particle with negative electric charge —Cj 
whose position at instant t is denoted by rj, this force is given, in gaussian units, by 



where E and B are, respectively, the electric and magnetic fields of the electromag- 
netic wave, c is the velocity of light in vacuum and v 3 - = ij is the velocity of the 
j-th particle at instant t. Denoting by Xj the difference between rj and the equilib- 
rium position of the particle, denoted by rj , we are lead to the following equation of 
motion for such a particle 





Xj + 2 7i Xj + ufe = 3 - Ej(t) + ^-x Bj(t) 



(2) 
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where, in order to simplify the notation, Ej(t) and Bj(i) stand for the electric and 
magnetic fields at instant t and position rj , that is, Ej(t) = E(rj ,t) and Bj(i) = 
B(rj ,t), and the quantities jj and ujj are the damping constant and the natural 
oscillation frequency associated to the j-th charged particle. In (J2]), we have assumed 
that |xj| is small enough so that the spatial variations of the electromagnetic fields 
over a distance of the order of the atom diameter can be neglected (|xj| <C A, with 
A being the wavelength of the incident electromagnetic wave). In other words, the 
electromagnetic fields can be considered as if they were uniform in space for each 
charged partice (this is the so called dipole approximation). In fact, since we will 
need to talk about the polarization P and magnetization M of the material, which 
are macroscopic quantities defined as volumetric densities computed over a volume 
5V which is very small macroscopically but big enough to contain many thousands 
of atoms, we shall assume that the condition |xj| 3 <C SVj <C A 3 is valid. 

For simplicity, we consider a linearly polarized and monochromatic plane wave of 
angular frequency uj and choose the cartesian axis in such a way that the electromag- 
netic fields are written as 

Ej(t) = E u costW — aij) y ; Bj(£) = E u cos(ut — acj) z ; k = kSt , (3) 

where aij is a constant, k = uj/c, x, y and z constitute a right-handed orthonormal 
set and E^ is the amplitude of the electromagnetic wave at position r j0 . It depends 
on Tjo, but this fact was not indicated explicitly to mantain the notation as simple 
as possible. Since the time averages to be taken do not depend on ay, we shall take 
otj = without loss of generality (alternatively, we can get rid off ay with a simple 
time translation t — > t + otj/oS). With this in mind, equation ([2]) takes the form (in 
components) 

Xj + 2'jjXj + uJjXj = — -— cos cut (4) 



c m 3 



Xj \ &iE(j_ 



Vj + 2 Wj + UJ jVj = ~ ( 1 — ) — - cos ut ( 5 ) 

Zj + 2'jjZj + ujjZj = (6) 

where Xj = ayx + yjy + zji. Looking for stationary solutions only, we see that the 
motion of the particle is restricted to the OXy plane, perpendicular to the magnetic 
field. Restricting to non-relativistic cases, that is, assuming \ij\ <C c and \ijj\ <C c, 
equations (j4]) and (JEJ) can be approximated, respectively, by 

p-E 

Xj + 2jjXj + ujjXj = and ijj + l^jijj + cu^yj = — - — — cos cut . (7) 

Hence, at this order, the stationary solutions of equations (@J, jSJ and are 

Xj[t) = ; yj (t) = cos (oot + <f> u ) ; Zj (t) = , (8) 



in which A w = Sj/Z^ and 



c _ e jE^ rv 
£"i — i A, 



2 2\ 2 , a 2 2 

UJ - - UJ ) + 



, cos <\) w = , sm <p u = — — . 

(9) 
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Observe that the particle oscillates in the same direction of the electric field with the 
same angular frequency u, but shifted by a phase 

Suppose that in a region with a volume SV (where the field variation is negligible 
by assumption) there are rij = NjSV oscillators with charge — ej, mass rrij and angular 
frequency Uj. Therefore, the total number of oscillators in SV is then N$v = J2j n j = 
J2j Nj5V and the polarization due to the present field is 

p = -Ete- ( 10 ) 

j 

The mean energy per volume in a period, reversibly stored in the medium, is 

U P = -<P • E) = E u £ ^(cos (ut + 4> u ) cosut) = - l -El £ M (^-) . 

j * j rrij \ ) 

(11) 

On the other hand, this same energy may be written as 

V P = -\p a E a = -\ X »K^ (12) 

where we used that P w = x' ( w ) E u for the Fourier transform of the polarization, with 
x'(u;) being the electric susceptibility of the medium. Comparing equations (ITT]) and 
(112]) . we identify 

x» = S , \Y — L i (13) 

where a, = e^/mjU) 2 is the static polarizability of the j-th oscillator. Comparing last 
equation to the following relation between polariazation and electric field in a linear 
and isotropic material (see, for instance, Chapter 4 of Ref. [I]) 

we obtain 

£ <^f^M => *0 = ! + d5) 

where e(u;) is the dielectric "constant" of the material (it is not a constant since it is 
a function of u). Last approximation, x' ( u ) ^ 1> i s valid specially for low density 
materials, as gases, for instance. Considering non-magnetic materials and neglecting 
the magnetic contribution, we have fi(uj) = 1, which makes the refractive index n(uS) 
approximately equal to 



n(u) = v / ^M~l + 27rx / H . (16) 



The behavior of n(u) = ye(u) is shown in Figure UJ for a single natural frequency 
Uq dominant in the medium. Between uj ~ ojq — 70 and u ~ u + 7o we see a region 
of anomalous dispersion (dn/du < 0). 
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The mean power per volume absorbed by the medium is given by 

N-e 2 N-e 2 -v-^; 2 

V{u) = - £ iV,e,(x,.E) = - E ZT^ ^(sm (ut + 4> u ) cos^t) = £ -JLL 2^- i£ 

(17) 

In terms of the mean value of the radiation intensity I u , the above equation takes the 
form 

reflecting the intuitive fact that the power is proportional to the intensity. 

In Figure [1] we show the behavior of ffl8l) as a function of the angular frequency 
with a fixed intensity. We considered also that the medium has only one dominating 
natural frequency u . We see that V{ui) is very small at all frequencies except those in 
the anomalous- dispersion region. Specially, for u = u>o the electric field is in resonance 
with natural oscillators of the medium and the power absorbed is maximum. 




Figure 1: Refractive index n{uj) = Je(u) (solid line) and absorbed power V(u) 
(dashed line), for fixed intensity I w as functions of the field frequency, for a medium 
with only one natural frequency u> (for convenience, we used the ratio uj/ujq). In 
both graphs we used a natural width 7 = 0,lu . The values on the vertical axis are 
written in arbitrary units. 



3 Magnetic contribution to dispersion 

In this section, we shall repeat the previous calculations but without neglecting the 
influence of the magnetic field. However, we shall compute only the first order cor- 
rection (in order |xj|/c) to the stationary solutions written in (jSJ). As it will become 
evident in a moment, to obtain the new stationary corrections to these equations we 
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need to obtain only the new stationary solution for the component Xj(t), namely, 
to consider the previously neglected iij/c term on the rhs of equation (j3J). Hence, it 
suffices to substitute yj(t) = A w cos(W + 4> w ) into equation (jll), which leads to 



Xj + 2'jjXj + uijXj 



-uAu sin(o;t + <p u 



rrij c 



cos(ut) 



~— sin(co>t + (j) u ) cos(ut) 
cZ^ 



2cZ u 
c Zl 



sin(2cj)f: + cj) u ) + sin <fr 
sin (2ut + 



2c Z, 



(19) 



where we used the definitions A u = Sj/Z^ = —ejE^/^rrijZ^), the trigonometric 
identity sin a cos b = (l/2)[sin(a + b) + sin(a — b)] and that sin^ = —l^jijjjZ^. 

A simple inspection on the rhs of the previous differential equation allows us 
to state that its stationary solution will oscillate with angular frequency 2a>, rather 
than lu. Following the usual procedure of obtaining stationary solutions of damped 
harmonic oscillators under an external force given by a constant term plus a harmonic 
one, we can obtain in a straightforward way the expression for the stationary solution 
for Xj(t). For convenience, we write below also the stationary solution for yj(t) up to 
first order in |x,- 



'c: 



Xj{t) = ^ukA 2 UJ -lkA UJ A 2uJ sm{2ujt + (j) ilJ + 



A w cos (cut + (j) u 



(20) 
(21) 



where k = u/c and the quantities A 2uJ , Z 2 ^ and 02^ can be obtained from equation 
dHJ) by just substituting u by 2u. As we shall see, equation (|2~T|) is, indeed, the correct 
expression up to this order. But before justifying this statement, a few comments are 
in order: 

(i) the stationary solution for zj(t) remains Zj(t) = 0, so that, as mentioned before, the 
stationary motion occurs on the OXy plane, which is perpendicular to the magnetic 
field; 

(ii) while the motion of the charged particle along the Oy direction is an oscillation 
around an equilibrium position at yj = 0, its motion along the OX direction is an 

oscillation around the equilibrium position Xj = -\ ukA^, which depends on the 

u 'j 

frequency of the electromagnetic fields. This dependence is similar as that of the 
absorbed power by the medium, namely: it is zero for uj = or uj — > oo and has its 



CO 



31 



maximum at uj 

(in) comparing the amplitudes of oscillation of equations (!20l) and ( |2T|) we see that 
the former is kA 2U ) times the latter. Hence, the oscillations in the direction of the 
electric field (Oy direction) are greater than those along the direction of the wave 
vector k (OX direction). 

Let us now estimate the error committed by neglecting the term Xj/c in (jSJ). Sub- 
stituting (120]) into (jSJ), and writing the solution of the resulting differential equation 
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as yj(t) = Vj(t) + Ayj(t), where y*j (t) is given by expression ([8]), it can be shown that 
Ayjit) is given by 



Ayj(t) = --k 2 AlA 2ul cos (ut + 2^ + (p 2oj ) - -k 2 A u] A 2u) A ZuJ cos (2ut + <\) w + 2w + 

(22) 

where the quantities A 3uJ and (p 3uJ can be determined in the same way as A 2u} and <p 2uJ , 
described before. Looking at the last equation, and remembering that k 2 = w 2 /c 2 , 
it is not difficult to conclude that Ayj(t) is of order |xj| 2 /c 2 and hence it can be 
neglected if we maintain terms only ut to first order in |xj|/c. That is why Ayj(t) is 
not present in equation (|2~T1) . 

The magnetic dipole momentum of the j-th oscillator is given by 



rri, 



Tc {Xj x Xj ' 



— CjSjk A U A 2(J] z 



cos(o;t + 4>w) cos(2c<;t + (f)^ + 4> 2u ) + 



+ cos{ujt + (p 2u) ) ^ — — sm(u;t + cp u 



(23) 



so that the corresponding medium magnetization is 



(24) 



and the medium mean energy per volume stored in the medium reads 
U m = -(M-B) 



N-e 4 E 4 



UJ 



, 4m3c 2 Z 2 Z,„. 



3 IjUJ Z 2uJ . 
-cos 02^ 5~— sin. 

4 u) i Z, 



E 



4m? c 2 Z 2 



3 /cj 2 -4cj 2n 
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7 2 



2^ 



2 V^Z W/ 



(25) 



Following a procedure totally analogous to that employed in Section 2 for the electric 
polarization, it is possible to write the magnetic energy as 



U„ 



ill 



V 8m?c 2 Z 2 



u 2 - ALU 1 



7 2 



2^ 



2 (l^u 

3 Vw i ^ / 



(26) 



where M = M u z is the mean induced magnetization parallel to the magnetic field. 
Comparing last equation to the relation between magnetization and the magnetic 
induction B in a linear medium, namely, (see, for instance, Chapter 5 of Ref. [4]) 



MM = — 



1 



4tt \h{u) + 2 



B, 



(27) 
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we obtain, after some algebraic rearrangements, 
3 + 8vrx ' ' 



3 - 4vrxr 



fj,(u) ~ 1 + Anxm (w) 





-4lu 2 


8 7 > 




[Jj - 4u, 2 ) 


2 + " 

+ 16 7 2 u; 2 


" 





(2f 



(29) 



In the last equation, we defined /3j = e|/m|c 3 u;| = a 2 /rrijC 3 and one could, in prin- 
ciple, think of Xm(w) as the magnetic susceptibility of the medium to the magnetic 
field, something analogous to the electric susceptibility x'( u ) defined by equation 
f|T2|) . Note that Xm(w) has a dependence on the frequency of the electromagnetic 
wave which is much more complicated than that exhibited by the electric susceptibil- 
ity x'( w )) as it is illustrated in Figure [H 




Figure 2: Eelectric and magnetic susceptibilities, given by equations (Tl2l and (1291 
as functions of the ratio uj/ujq. The dashed line stands for while the solid line 

corresponds to Xm{w)- F° r simplicity, we considered only one natural frequency, Uq, 
(and the corresponding natural line width 70). We have made 70/^0 = 0.1 in both 
graphs. The values on the vertical axis are written in arbitrary units. 



Differently from x'( u )i which varies smoothly with the frequency except in the 
range of anomalous dispersion (with a width 270 centered at uq), x m (w) shows an 
anomalous behavior around u = uj /2. There is also a pronounced peak around 
uj = ujq, similar to that appearing in the absorbed power by the medium, as shown in 
Figure [H but whose intensity varies with I/7 2 instead of l/ 7 o- 

In Figure [3] we plot some graphs for ^ m (w) versus uj/ujq for different values of 7 . 
In the range \ojq < u < u we observe an uncommon behavior with the frequency, 
while for high frequencies, uj 3> oj , ^ m (w) goes to zero faster than x'{ u )- It is worth 
mentioning that, for low frequencies, uj <C uj q , only x m ( w ) goes to zero. 

In contrast to the electric permittivity e(u>), given by equation (Till , the magnetic 
permeability fi(uj) depends on the intensity I w of the electromagnetic wave, as can 
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Magnetic Susceptibility vs y 




Figure 3: Magnetic susceptibility x m ( w ) as a function of the ratio uj /ujq for different 
values of 70 . The solid line corresponds to 7o/cu = 0.08, while the dashed line, to 
7o/co> = 0.1, and the dotted line, to 7o/^ = 0.12. The values on the vertical axis are 
in arbitrary units. 



be seen from equations (|28|) and (|29|) . For x m (w) <C 1, which is generally the case, 



(j,(oj) — 1 is proportional to I u 
index takes the form 



Including the magnetic contribution, the refractive 



n{uj) 



n(uj) ~ 1 + 2rr 



(30) 



We see that the role of the magnetic field is to sum 2iiXm(^) to the rhs of equa- 
tion (TL6I) . which introduces a linear dependence with the intensity of the electromag- 
netic wave for fixed frequencies. This term, however, is usually much smaller than 
the electric contribution, but may become relevant for values of intensity that make 
X'(uj)/xm(u) ~ 1- in the next section we discuss if there is a possibility of check- 
ing experimentally this magnetic contribution. Figure H] shows the effect of this new 
magnetic term on the expression for the refractive index. 

For low values of the intensity I w , but still high enough to make the contribution 
due to XmM a relevant one, the frequency dependence of the refractive index is very 
similar to that shown in Figured], except for the appearence of a range of a secondary 
anomalous dispersion (SAD), centered at uj = ^uj (besides the usual anomalous 
dispersion around uj = uj ). This fact may be viewed as a first signature of the 
influence of the intensity on the refractive index. 

As the intensity increases, the behavior of n(uj), in the interval ^uj < oj < uj , 
becomes more and more similar to that of x m (c<j). We also notice that the intersection 
between the curve for n{uj) and the line n(uj) = 1, which occurs at uj = ujq for vanishing 
intensity, is shifted as the intensity is increased. Finally, for uj <ti uj and uj > uj , 
n(uj) does not change appreciably when one varies the intensity. 

In Figure [5] we see the relation between n{u) and frequency for different values 
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Refractive Index vs Intensity 




Figure 4: Refractive index, given by equation fl30|) . for different intensities of the inci- 
dent electromagnetic wave and 70/^0 = 0.1o>o- The curve with solid line corresponds 
to an intensity 5 times smaller than the intensity used to draw the curve with dashed 
line and 10 times smaller than that used in the curve with dotted line. The values on 
the vertical axis are written in arbitrary units. 



Refractive Index vs y 
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Figure 5: Refractive index, given by ( 1301) . as a function of uj/ojq for different values 
of 70 and fixed intensity. The solid line stands for 70/cuo = 0.03, the dashed line, to 
7o/cg> = 0.05 and the dotted line, to To/^o = 0.1. The values on the vertical axis are 
in arbitrary units. 



of To, keeping the incident intensity fixed. As expected, as To becomes smaller and 
smaller, the secondary anomalous dispersion becomes more and more evident in the 
interval u = ^cuq. 



10 



4 Numerical estimatives 



So far we have discussed the problem of the interaction between microscopic struc- 
tures and electromagnetic radiation until first order in v/c, without worrying about 
numerical values that could sustain the possibility of experimental observation of 
those predictions (magnetic corrections to the reactive properties of substances). In 
this section, we want to analyze the limits for the parameters associated to the ra- 
diation (frequency and intensity) and to the system (natural frequencies) that make 
such predictions observable within the validity domain of this model. 

Henceforth, we consider our material oscillators as electrons with electric charge 
e = — 4,803- 10~ 10 (erg • cm) 1//2 and mass m = 9,109- 10 _28 g. We shall call them actives 
since they are the only charges that can interact with the incident radiation. For the 
sake of simplicity, we shall assume that there is only one natural frequency present, 
denoted by ujq, so every electron oscillates with the same frequency. Therefore, the 
sum in (1121) . (118j) and (1291) keeps only one term. We also admit that these oscillators 
are perfect harmonic oscillators, an assumption that leads to 

* S 5S?" (31) 

Strictly speaking, this model is applicable only to substances formed by active elec- 
trons that behave as harmonic oscillators. However, it can be successfully extended 
to those that actually are not formed by oscillators, but in these cases the parameters 
atj, j3j, ojj and jj have to be considered independent from each other, being adjusted 
through experimental data. 

Let us, then, analyse the secondary anomalous dispersion. In the range where it 
occurs, we can approximate the magnetic susceptibility Xm(w) by 

3 [u - 2uj) + 7g 

while, in this range, the electric susceptibility is almost constant, and can be taken 
as being 

x'(w)~^W (33) 

We have seen that magnetic effects will become relevant for \x'(uj)/xm{^)\ ~ 1- The 
maximum value of |% m (u;)| occurs around = (c^o ± 70) and is given by 

77 N 8 T U ° - nN ° e2 
— iV po-'w — — — 5 — = 

3 70 m z ujQ 

so that the minimum intensity needed is 

fiV a 327r 2 mc 3 //im\ 3 2 

~ pww^ = I" IT K 260 (17) GW/cm ' (35) 

These intensities are relatively small in the microwave region while they are very high 
in the optical region: for Ao = 0,5/im the intensity is something like 2,0TW/cm 2 , 
making it hard to observe the secondary anomalous dispersion in that range of fre- 
quencies. Nevertheless, in the range from infrared to microwaves, an experimental 
observation of these magnetic effects seems feasible: for Ao = 100/im, for instance, 
260kW/cm 2 , a value that can be achieved nowadays by lasers. 



Xm(w) ^ — Nq(5qUJqI w — ° o , 2 , .2 ' ( 32 ) 



\Xm\ ^ o N 0PoIu— = ~ 2 5 ' ( 34 ) 



11 



5 Conclusions and final remarks 



In this work we have discussed the Lorentz model for dispersion and absorption of 
electromagnetic radiation in diluted, linear and isotropic material media. We started 
reviewing the usual computation of the electric permittivity of the material that 
describes the response of the material oscillators only to the electric field of a linearly 
polarized and monocromatic electromagnetic wave. In this simplified model, the 
magnetic field is not taken into account, since its contribution is already of order v/c. 
From this model, one obtains the main results about the electromagnetic reactive and 
dissipative properties of matter, namely: the dependence of the refractive index with 
the radiation frequency (dispersion); the existence of a region of strong anomalous 
dispersion for the frequency range \ujj — uj\ ~ 7j (which we call, for convenience, a 
primary anomalous dispersion), where dn/du < 0; the absorbed/dissipated power is 
proportional to the intensity of radiation and varies quickly in the range | c^j — | ~ jj, 
achieving its maximum value at resonance. 

Then we discussed what are the first corrections to the previous model when we 
include the magnetic term of the Lorentz force into the equations of motion of the 
material oscillators. Despite this term does not change substantially the dissipative 
properties of the medium, it leads, when solved up to first order approximation in 
v/c, to quite interesting (and presumably observable) effects concerning the reac- 
tive properties of matter, namely: (i) the refractive index acquires a dependence 
on the intensity of the electromagnetic wave (besides its dependence on the wave 
frequency); (ii) the appearance of a secondary zone of anomalous dispersion in the 
region — u\ ~ jj and of an intensification of the primary dispersion peak. This 
result breaks the monotonic increasing behavior of the refractive index for sufficiently 
high intensities, in the range \ojj < uj < ujj] 

The observation of the new magnetic effects just described in the resonance region 
(uj ~ ujj) is frustrated due to the strong energy absorption. However, the situation 
may not be the same in the secondary anomalous dispersion zone, where absorption 
is very small, so that the required high intensities of the incident radiation may be 
achieved. In order to investigate the plausibility of making real experiments with the 
present technology with the purpose of observing the secondary anomalous dispersion, 
we have presented in the previous section some numerical estimatives. The observa- 
tion of this secondary anomalous dispersion would be a signature of the dependence 
of the refractive index on the intensity of the electromagnetic wave. 

As a final comment, we would like to say that, so far, we have considered only 
a linearly polarized incident radiation. However, one can use circularly polarized 
radiation as well. In this case, though the magnetization acquires a non-vanishing 
component along the direction of propagation of the electromagnetic wave, the final 
results remain the same, since the direction of propagation is perpendicular to the 
magnetic field, so that the contribution of this extra term vanish in the expression of 
(M-B). 
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